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phase separated state of Lai ^Ca^MnOs 
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The chemical influence in the phase separation phenomenon that occurs in perovskite manganites 
is discussed by means of ab initio calculations. Supercells have been used to simulate a phase 
separated state, that occurs at Ca concentrations close to the localized to itinerant crossover. We 
have first considered a model with two types of magnetic ordering coexisting within the same 
compound. This is not stable. However, a non-isotropic distribution of chemical dopants is found to 
be the ground state. This leads to regions in the system with different effective concentrations, that 
would always accompany the magnetic phase separation at the same nanometric scale, with hole-rich 
regions being more ferromagnetic in character and hole-poor regions being in the antiferromagnetic 
region of the phase diagram, as long as the system is close to a phase crossover. 

PACS numbers: 



I. INTRODUCTION 

Perovskite manganites belong to a class of mate- 
rials that have become increasingly important in the 
last few years, particularly because of the large va- 
riety of phenomena they exhibit, closely tied to the 
strongly correlated nature of their electronic struc- 
ture. The series of manganites Ai^B^MnOs (A 3+ 
being a lanthanide and B 2+ an alkaline earth) is 
one of the most intensively studied families X The 
interest in these materials grew significantly after 
the discovery of colossal magnetoresistance (CMR), 
a large negative magneto-resistance taking place un- 
der an applied magnetic field; 2 - but also considerable 
research has focused on their unusual properties such 
as phase separation, 3 charge and orbital ordering 4 ^ 
and also for their potential applications in thin-film 
devices^— derived of their peculiar electronic struc- 
ture properties. 

The case of A = La 3+ and B = Ca 2+ has become 
a prototypical example of manganite compounds 
showing a very rich phase diagram.! 11 : 12 The region 
in the phase diagram with concentrations of about 
18% and 50% Ca 2+ have been thoroughly stud- 
ied, because samples exhibit phase separation, in- 
trinsic to the proximity to the magnetic/electronic 
phase transition. In that part of the phase diagram, 
Lai-^Ca^MnOa presents mixed Mn 3+ /Mn 4+ valen- 
cies and a strong electron-phonon coupling, the mag- 
netic coupling being a modified version of the dou- 
ble exchange interaction to take into account the 
coupling between the electronic states and vibra- 



tional modes (vibrons). The mechanism for ferro- 
magnetic (FM) ordering in manganites exhibiting 
CMR has been recently considered^ suggesting a 
breakdown of the adiabatic approximation for the 
itinerant electrons. These results are related with 
dynamic phase segregation into hole-rich itinerant- 
electron regions and hole-poor localized-electron re- 
gions, leading to suppression of the conventional 
double-exchange mechanism. 

Several theoretical studies have been performed 
in these manganite perovskite systems using dif- 
ferent techniques. The electronic structure of 
Lai-ajCa^MnOa was calculated using the LAPW 
(Linearized Augmented Planewave) methodi 14 : 15 
Transport properties were studied using the virtual 
crystal approximation, neglecting disorder effects in 
the latticed The formation of a pseudogap follows 
from numerical calculations based on the dynam- 
ical mean field theory, indicating the existence of 
a charge-nonuniform ground stated An extensive 
study of the physics of manganites, including the 
structure we are studying here, was carried out by 
Salomon et ali^ 8 - Unfortunately, an ab initio descrip- 
tion of the phase separation in manganites has not 
been given in the literature. 

It is now widely accepted that some sort of elec- 
tronic phase segregation is at the origin of the CMR 
effect in manganites. Moreo et aZ.— proposed a sce- 
nario where the coexistence of two phases with differ- 
ent charge density is unstable against breakup into 
nanometer size clusters due to the Coulomb force. 
Two different kinds of electronic phase segregation 



have been proposed to exist in manganites:™ equal 
density vs. different density electronic phase seg- 
regation at the micro- and nano-scale, respectively. 
Experimentally, magnetic phase separation has been 
observed in the form of nanoscale phases with differ- 
ent magnetic properties^ i.e. regions in the crystal 
with nanometric size that present a different mag- 
netic ordering. Imaging techniques 2 ^ have recently 
shown the close relationship between nanoscale 
phase separated regions of an electronic/magnetic 
origin and CMR. According to those studies, these 
regions, on the several nm range, have not chemical 
origin. A different view comes from an analysis of 
the influence of the dopant distribution in the super- 
conductivity and phase separation that has recently 
appeared 2 ^ studying the chemical inhomogencitics 
at the nanoscale in (Bai_, c K. I ;)Fe2As2 superconduct- 
ing pnictides. 

The fundamental problem is to give a correct de- 
scription of how the holes are distributed in the sys- 
tem. It has been known for long that even the sim- 
plest models 2 ^ predict a non-uniform density distri- 
bution of dopants, which are unstable against seg- 
regation into hole-rich and hole-poor regions. To 
gain insight into this problem, in this paper, we 
have taken into account a structure with Ca 2+ dop- 
ing levels close to x ~ 0.2 (exactly x = 0.1875) 
of Lai-zCa^MnOs from a computational point of 
view via "ab initio" calculations. Typically, this is 
not done, and models are utilized instead, due to 
the large supercells that are required, but even the 
crudest supercells can help us, as we will see, to un- 
derstand the electronic and chemical origins of phase 
separation. It is worth recalling that the above men- 
tioned doping level is experimentally placed around 
the critical boundary in which the FM insulating to 
FM metallic transition occurs at low temperatures 
with respect to doping . 11 ' 12 We have constructed 
several superstructures in different magnetic config- 
urations to explore the magnetic origin of phase sep- 
aration. Since our calculations yield a homogeneous 
FM ground state, we will analyze the influence of 
the Car + dopant distribution in the lattice in this 
magnetic ground state. Varying the Ca 2+ positions 
from a random distribution to a slightly correlated 
one, we can analyze the influence of dopant distri- 
bution in the electronic structure and the physical 
properties of the compound by means of "ab initio" 
calculations. 

For the sake of completeness, we will make a sim- 
ilar analysis in a very different Ca 2+ concentration 
(x = 0.375), in which the ground state is FM metal- 
lic becoming paramagnetic insulating at high tem- 
peratures. The idea of that calculation is to study 
a concentration away from the phase boundaries at 
x = 0.1875 and also from x — 0.5, where the ground 



state becomes an insulating antiferromagnetic (AF) 
and charge ordered (CO) state. At those two hole 
densities there is a strong tendency towards phase 
separation caused by intrinsic inhomogeneitiesi but 
x = 0.375 (unlike x = 0.1875 and 0.5) is not close 
to the boundary between different ground states in 
the Lai-zCa^MnOs solid solutionj 11 ' 12 These anal- 
yses will help us discuss the consequences of having 
chemical inhomogeneities on the nanoscale in these 
phase separated compounds. 



II. COMPUTATIONAL DETAILS 

Electronic structure calculations were performed 
within the density functional theory framework 2 ^ us- 
ing WIEN2K software ; 26 ' 27 which utilizes an aug- 
mented plane wave plus local orbitals (APW+lo) 28 
method to solve the Kohn-Sham equations. This 
method uses an all-electron, full-potential scheme 
that makes no shape approximation to the poten- 
tial or the electron density. 

For modeling the correlated behavior of the d elec- 
trons of the system, we have included strong cor- 
relation effects by means of the LDA+U scheme^ 
controlled by an effective U (U e //=U-J), being U 
the on-site Coulomb repulsion and J the on-site ex- 
change constant (taken as J=0, as is common prac- 
tice in literature). 30 Values of U e //: 2.7, 4.0 and 
5.5 eV were used in the calculations. These values 
are on the range that have been used in the past 
for various Mn oxides . 31 ' 32 The structural minimiza- 
tion was carried out using the GGA— PBE scheme 3 ^ 
(generalized gradient approximation in the Perdew- 
Burke-Ernzcrhof scheme) minimizing the forces on 
the atoms and the total energy of the system. The 
parameters of our calculations^ depend on the type 
of calculation but for any of them we converged with 
respect to the fc-mesh and to RmtKmax, up to 75 k 
points (15 in the irreducible Brillouin zone, in a 3 x 
3x5 sampling of the full Brillouin zone for a su- 
percell N x N x N perovskite unit cells) and up to 
RmtK max = 6.0. Muffin tin radii chosen were the 
following: 2.43 a.u. for La, 2.29 a.u. for Ca, 1.96 
a.u. for Mn and 1.73 a.u. for O. 



III. RESULTS AND DISCUSSION 

A. Magnetic study 

The goal of the first part of our work is to simulate 
a magnetically phase separated state by embedding 
a FM (AF) phase into an AF (FM) matrix. For this 
reason, we have set up several superstructures based 
on the unit cell Lao.8i25Cao.i875Mn03 (at a concen- 
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FIG. 1: (Color online) Schematic cases of the phase sep- 
arated structures studied. On the bottom left we see a 
representation of the unit cell of Lao.8i25Cao.i87sMn03 
(lxlxl) in the ferromagnetic (FM, blue) configuration 
and in the antiferromagnetic (AF, red) configuration (a). 
On the upper left (b) we represent the (2x2x1) cases, 
one dimensional AF chains surrounded by a FM guide 
and vice versa. On the right (c) we represent both realis- 
tic cases of phase separation (2x2x2): an AF nanomet- 
ric region completely surrounded by a FM region (upper) 
and vice versa (lower). In this latter case, a 4x4x2 is 
depicted for clarity. 



tration close to a metal-insulator transition in the 
phase diagram), using a 2x1x2 and a 2x2x2 su- 
percell (in perovskite unit cell units), as shown in 
Fig. [TJ The larger and initially more realistic cases 
(2x2x2) in Fig. [T] (c), serve as a model of the ex- 
perimentally observed situation of a magnetic phase 
(antiferromagnetic (AFREAL) in the upper figure, 
ferromagnetic (FMREAL) in the lower one) com- 
pletely surrounded by a magnetically different phase 
(FM and AF, respectively). Figure [1] (b) shows the 
cases of an AF chain surrounded by a FM guide 
(AFCHAIN) and vice versa (FMCHAIN). We have 
also calculated the system in the simple ferromag- 
netic (FM) and G-type AF cases (Fig. 00(a)). 

We have calculated all these phase separated su- 
perstructures and compared their relative stability 
(total energy calculations) with respect to an en- 
tirely FM structure and also an entirely AF struc- 
ture (G-type). Results are summarized in Fig. [2j 
We can see that the purely FM structure is the most 
stable one in every case, the stabilization energy 
ranging from 130 meV/Mn compared to the FM- 
CHAIN and AFREAL structures to 400 meV/Mn 
with respect to the AF structure. Calculations pre- 
dict that these magnetically phase separated states 
are not stable. The picture is the same for all val- 
ues of U we have calculated. From these calcula- 
tions, we can conclude that the FM configuration 



FIG. 2: (Color online) Calculated energies of the con- 
structed structures at several values of U within the 
LDA+U approximation. We take the more stable cal- 
culated case (FM case) as the zero energy state. All 
energies are expressed in meV/Mn. 



at the concentration xq — 0.1875 is the most stable, 
in agreement with experimental observations at that 
concentration ! 1 1 i 12 These results suggest that the co- 
existence of two magnetic phases in the compound 
require some other changes in the crystallographic 
and electronic structure for it to be stable. 



B. Chemical study 

Holes in the e g Mn levels can be introduced by the 
substitution of La 3+ by Ca 2+ over the entire com- 
positional range in the Lai-^Ca^MnOs solid solu- 
tion. Additionally, hole doping can be created by 
means of cationic vacancies. Lao.si25Cao.i875Mn03 
is a distorted perovskite with a pseudocubic lat- 
tice parameter (a c ) of 3.89 A. The tilting of the 
MnOe octahedra results in an orthorhombic symme- 
try (space group Pnma, no. 62) with lattice param- 
eters a w c sa y/2a c and b a s» 2a c for this composi- 
tion. The ability to oxidize Lai-zCa^MnOa materi- 
als with x < 0.2 leads to important structural, elec- 
tronic and magnetic variations related to the pres- 
ence of cationic vacancies . 10 ' 35 ! 36 Thus, differences 
between non-stoichiometric and stoichiometric sam- 
ples should be always considered. 

Stoichiometric Lao.8i2sCao.i875Mn03 compound, 
i.e., with no cationic (and anionic) vacancies is con- 
sidered in the present study. We initially have as- 
sumed that the La 3+ and Ca 2+ atoms (the dopants) 
are distributed in a perfectly homogeneous manner 
throughout the crystal. But, experimentally, this 
might not be the case; nano-sized chemical inhomo- 
geneities in the structure can happen and will not be 
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FIG. 3: (Color online) Schematic representa- 
tion of Ca 2+ atoms distribution in the lattice of 
Lao.8i25Cao.i875Mn03. <r> is the average distance 
between Ca 2+ atoms in the lattice. These figures show 
the four possible atomic ordering of Ca 2+ atoms in the 
lattice (because of its symmetry) and its <r> from 
the most inhomogeneous case distribution (a) to the 
homogeneous case distribution (d). 

detected below the few nm limit (22 but they could 
have a key role in the properties of a phase separated 
state. 

The next step of the work is to study if such 
dopant inhomogeneities can occur in the nanoscale. 
To do this, we have varied the atomic positions of 
the Ca 2+ atoms to distribute them in various man- 
ners throughout the crystallographic unit cell. In 
the Lao.8i25Cao.i875Mn03 structure we can change 
the distribution of the Ca 2+ cations to create four 
possible inequivalent chemical configurations. These 
can be seen in Fig. |3j We can define the quantity 
<r> as the average distance between Ca 2+ atoms in 
the lattice, and use it as a measure of how homoge- 
neously these dopants are distributed: 

»=i 

where di is the distance between nearest Ca 2+ 
atoms and N is the number of nearest Ca 2+ atoms. 

In principle, inserting chemical inhomogeneities 
into the lattice will generate uncompensated forces 
and increase the elastic energy in the structure, but 
we need to calculate all the energetic terms involved 
to obtain the most stable configuration. In our case, 
we will do that using ab initio calculations. 

We have constructed these four possible cases us- 
ing a FM ordering, the most stable magnetic struc- 
ture we explored previously. We have performed a 



<r> 

FIG. 4: (Color online) Energies of the calculated struc- 
tures for Lao. si25 Cao.i875Mn03 (y axis) in the FM con- 
figuration. The energies are normalized to the lowest 
value of the four possible Ca 2+ chemical distribution con- 
figurations. On the x axis we write the values of <r> 
for each structure (see text), from the larger one on the 
right (total homogeneity structure corresponds with Fig. 
[3]d) to the shorter one on the left (total inhomogeneity 
structure corresponds with Fig. [3] a). The inset shows 
an analogous study of the second magnetic stable con- 
figuration calculated in section A of this paper, which is 
the FMCHAIN configuration. 



full lattice relaxation for each of the cases, and then 
we have calculated their total energies. The results 
are presented in Fig. [4] We can observe that the 
most stable structure corresponds to an intermedi- 
ate case between the most homogeneous (Fig. [3ji) 
and the most inhomogeneous (Fig. [3k.) case, which is 
the one used for the calculations in previous sections. 
Taking the ground state structure as the <r>= 6.71 
A, we found this one is more stable by about 19 
meV/Mn with respect to the second most stable 
configuration, and by 219 meV/Mn to the most un- 
stable configuration (largest chemical homogeneity 
attainable at this concentration) . This means that a 
non-homogeneous dopant distribution occurs in the 
material that would lead to a nanoscale distribu- 
tion of effective concentrations. Analogous results in 
terms of dopant distribution energetics can be found 
in other magnetic configurations (hence, it is not an 
issue of the FM solution). 

The most stable < r > value is larger than one 
cubic perovskite lattice parameter. Several per- 
ovskite subcells will be involved in forming this 
ground state Ca 2+ distribution. Moreover, varia- 
tions of the atomic position of the various Ca 2+ 
cations will imply different environments for each 
of the Mn atoms in the unit cell. Thus, the ex- 
istence of areas with such a typical size suggests 
the magnetic and electronic characteristics differ in- 
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FIG. 5: (Color online) Majority spin (positive values) 
and minority spin (negative values) partial density of 
states (DOS) plots of the eight inequivalent Mn atoms 
in the ground state (both chemically and magnetically) 
configuration. Fermi energy is represented by a vertical 
line at E=0 eV. 



side and outside the clusters leading to an electron- 
ically inhomogeneous state on a nanometric scale in 
Lao.8i25Cao.i875Mn03. Zones where Ca 2+ atoms are 
grouped lead to a local concentration displacement 
towards the hole rich region (Ca 2+ -rich) in the phase 
diagram of La 1 _. Z ;Ca ;r Mn03. Instead, zones where 
the density of Ca 2+ atoms are lower, we have a La 3+ - 
rich (electron-rich) zone, being a local concentration 
displacement to the AF region of the phase dia- 
gram. This scenario is consistent with the stable FM 
order' 55 near Mn 4+ , localized near Ca 2+ -rich regions, 
in agreement with the hole distribution described by 
Alonso et at as the hole attractor model. 36 Our 
calculations show that the existing chemical inho- 
mogeneities can be a coexisting factor that works 
together with the electronic/magnetic phase separa- 
tion. 

We have calculated the electronic structure of the 
material in the ground state (both in terms of mag- 
netic ordering and chemical configuration) . The val- 
ues obtained for the total magnetic moment in this 
structure are of 3.8125 [Ib per Mn atom, as ex- 
pected. This corresponds to an average Mn valence 
of +3.1875. We can also analyze the magnetic mo- 



ments inside the muffin-tin spheres for each of the 
Mn atoms. These are summarized in Table [H Their 
magnetic moments present differences on the order 
of 1 % at most. Also the total d charge inside the 
Mn muffin-tin spheres varies by only 0.015 electrons. 
This all would be translated in changes in effective 
Ca concentration x of only 0.01. This can be seen 
in the largely homogeneous electronic structure ob- 
served in the partial density of states of the inequiv- 
alent Mn atoms in the ground state of the compound 
(Fig. [5]). In a metallic solution, the way the differ- 
ent local concentrations correlate with an electronic 
phase separation in the system is not obvious, but 
we observe that the chemical inhomogeneities relate 
to a different electron count in each Mn atom. 



TABLE I: Values of the atomic magnetic moment in 
units of /is inside the muffin-tin sphere for the ground 
state configuration. Values of total magnetic moment 
are 3.8125 /is per Mn atom, with an average valence 
Mn 3 1875+ , the value expected from the electron count 
in a fully polarized half-metallic ferromagnet. 
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The typical sizes of these hole / electron rich zones 
that turn out from our calculations for the ground 
state chemical distribution can be seen in Fig. [6] 
We have calculated the effective x for each inequiva- 
lent Mn atom in the supercell according to its near- 
est neighbours (up to fourth neighbours), assum- 
ing the usual Ca 2+ /La 3+ valencies, for this most 
stable chemical configuration. We observe there a 
6.5 — 7.4 A diameter zone where the effective dop- 
ing of Lai_ x Ca x Mn03 has more Ca 2+ than the real 
macroscopic average concentration (x = 0.1875) and 
a 4.3 — 6.5 A diameter zone where the effective dop- 
ing is displaced to the La 3+ -rich zone in the phase 
diagram. These calculated cluster sizes are com- 
parable, but somewhat lower than those reported 
in literature 3 ^ concerning the Ca-rich region in the 
Lai_ x Ca x Mn03 solid solution. They show the im- 
portance of chemical ordering at the nanoscale when 
the system is close to a metal-insulator transition. 
The presence of chemical inhomogeneities on a sim- 
ilar nanometric scale than the experimentally ob- 
served electronic/magnetic phase separation could 
act as a precursor for the phase separated scenario to 
develop in the system. The largest deviation around 
the Xq =0.1875 in terms of Ca concentration (x) is 
about 0.04 (slightly larger than the equivalent devi- 
ation in the magnetic moments of Mn atoms, which 
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FIG. 6: (Color online) Difference between lo- 
cal concentrations and global concentration of 
Lao.si25Cao.i875Mn03. Orange (lighter) colors in- 
dicate a hole rich region on the phase diagram of 
Lai-zCa^MnOa,— leading to a concentration where 
FM is the ground state. Blue (darker) colors indicate 
an electron rich region on the phase diagram (where 
the ground state is AF). The diagram shows zones of 
about 7.0 A in typical size with Ca concentration x e ff 
> xo and other zones of about 5.5 A with x e ff < xo ( 
x = 0.1875 ). 



was calculated to be about 0.01). Chemical vari- 
ations are stronger than the electronic variations. 
Changes in the local doping of each Mn will vary 
the effective Ca concentration from 0.167 < x < 
0.208, which are enough to move the system to the 
other side of the magnetic/electronic phase transi- 
tion, that occurs at about x = 0.18. In fact, different 
nanosize domains have been experimentally detected 
around this compositional range: FM clusters of 16 
A for x — 0.17 whereas an approximate size of 8 A 
has been observed for x — O.2. 21 - This last result con- 
firms the existence of FM clusters confined in about 
two perovskite subcells, in a similar scale than the 
chemical inhomogeneities we find in our calculations 
for x — 0.1875. The zones with a lower Ca 2+ concen- 
tration will have an AF ground state and the zones 
with higher values of x will be in their FM ground 
state, according to the phase diagram of the system. 

In order to further check the reliability of these 
calculations, we have analyzed a different Ca 2+ 
concentration (x — 0.375) which is not at the 
boundary between different ground states in the 
Lai-zCa^MnOa solid solution. For this, we have 
constructed four different possible Ca-dopant distri- 
butions in the Lao.625Cao.375Mn0 3 i 38 i 39 A similar 
energy analysis as the one performed above shows 
that the chemical inhomogeneities that occur at this 
concentration are even larger in magnitude. The dis- 
placements in the local concentration of the < r > = 




<r> (A) 



FIG. 7: (Color online) Analogous study of energies for 
Lao.625Cao.375Mn03 than Fig. [4] Calculations were 
made in a FM configuration (the more stable at this 
concentration) at several values of U, leading to similar- 
results than Lao.8i25Cao.i875Mn03. 

5.5 A structure, that can be observed in Fig. [7j 
ranges from 0.335 < x < 0.407. Note that < r > 
is smaller than in the x = 0.1875 case (i.e. 6.71 
A ). This is related to the richer calcium concen- 
tration in x = 0.375. The above range of Ca 2+ 
concentrations is always far enough from the FM 
insulator ground state at x = 0.18 as well as the 
AFM with CO ground state observed at x = 0.5. 
In spite of this, experimental evidence of inhomo- 
geneities in manganites has been obtained. 1 Thus, 
chemical inhomogeneities are always present in this 
kind of samples, but a proximity to a phase bound- 
ary is required for them to act as nucleation centers 
of an electronic/magnetic phase separation scenario. 



IV. CONCLUSIONS 

Our results show important evidences about an 
often forgotten factor (mainly due to the compli- 
cations for its direct calculation and measurement) 
in studying the origin of electronic and magnetic 
phase separation close to a magnetic phase transi- 
tion. By means of ab initio calculations, we have 
studied mixed magnetic phases but none of them is 
stable. However, chemical inhomogeneities caused 
by cation disorder in a scale smaller than 1 nm are 
stable in the system according to our total energy 
calculations. Close to a FM/AF phase transition, 
these nanoscale doping inhomogeneities are large 
enough to produce hole- rich/hole-poor regions in the 
compound at the nanometer scale. For some con- 
centrations, like the x ~ 0.2 phase boundary, they 
are found to have a similar size than the observed 
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phase separated regions. Thus, their effects cannot 
be ignored and could be substantial, whether driving 
electronic/magnetic phase separation or accompany- 
ing it. Those cationic inhomogeneities could act as 
charge attractors and/or nucleation centers for the 
phenomenon of phase separation to develop. If the 
scale of the phase separation phenomenon were mi- 
crometric, or much larger than 1 nm, these chemical 
effects could possibly be ignored. 
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